Circuit-QED analogue of a single-atom injection maser: Easing, trapping states and 

multistability 
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We study a superconducting single-electron transistor (SSET) which is coupled to a LC-oscillator 
via the phase difference across one of the Josephson junctions. This leads to a strongly anharmonic 
coupling between the SSET and the oscillator. The coupling can oscillate with the number of photons 
which makes this system very similar to the single-atom injection maser. However, the advantage of 
a design based on superconducting circuits is the strong coupling and existence of standard methods 
to measure the radiation field in the oscillator. This makes it possible to study many effects that 
have been predicted for the single-atom injection maser in a circuit quantum electrodynamics setup. 



Introduction The ability to fabricate and control su- 
perconducting quantum circuits has given birth to the 
field of circuit quantum electrodynamics (circuit QED). 
Within these circuits, quantized charge, superconducting 
phase difference and even individual microwave photons 
can be controlled and manipulated, introducing the idea 
of 'artificial' atom-photon physics. Such circuits have 
been used to reach the strong coupling regimeii^, observe 
heating and cooling^Ti^, realize a three-level-laseri^— and 
make great strides forward in creating highly tunable 'ar- 
tificial' atoms coupled to an oscillator—. A number of 
proposalsiiii^ and ultimately experiments^!^ have real- 
ized the idea of a single 'artificial-atom' laser (or more 
precisely maser). Such work strives towards a circuit 
analogy of the experimental generation of coherent mi- 
crowave radiation using single-atom masers^^"— . 

Despite much progress, several important physical ef- 
fects that have been painstakingly studied in single-atom 
injection masers are yet to be realized in circuit-QED. 
These effects include multistability of the cavity fieldi^, 
trapped photon numbers states^i and the creation and 
annihilation of single photons^^. In all these experiments, 
the state of the system is inferred via measurement of the 
excited state occupancy of Rydberg atoms as they drop 
through a microwave cavity. 

In this paper we will investigate a superconduct- 
ing single-electron transistor (SSET) coupled to a LC- 
oscillator via the phase difference across the SSET's 
right-hand Josephson junction (JJ). In doing so, we show 
that this device allows one to reach a regime of strongly 
nonlinear coupling. It displays multistability and trap- 
ping states, as well as possessing an operating point 
which provides protection against low-frequency charge 
noise. As such, this device provides a close analogy to the 
single-atom injection experiments of cavity-QED. How- 
ever in contrast to previous experiments, a circuit-QED 
setup allows the direct detection of the photon state 
of the resonator via time resolved measurements of the 
emitted microwave radiatiorti^. 

The System In the following we discuss two possibe 
realizations of our desired circuit. One realization is the 
series coupling of the oscillator with the SSET, Fig. [T^. 
The other option is the more conventional inductive cou- 



pling scheme^. Fig. [T}d. In the latter it is necessary that 
a SQUID forms the effective right-hand JJ of the SSET. 
Although incorporating a SQUID will provide additional 
tunability to the former case as well. We divide the co- 
herent Hamiltonian into three parts, describing the artifi- 
cial atom, the photon mode and the interaction between 
them 
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In both realizations the artificial atom, that forms the 
basis of the micro-maser, is provided by the charge on 
the SSET island and Cooper-pair tunneling across the 
left J J, and is described by the Hamiltonian 



-ffatom = 4.EciN ~ Ncf - EjL cos(0l). 



(2) 



Here Ejl is the corresponding Josephson coupling and 
the charge on the island 2eN and the phase difference 
across the left J J satisfy the periodic commutation 



relation [N,e^'"^^] = ±( 
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The explicit forms of the 



gate charge Nq and the charging energy Eq depend on 
the realization and are given below. The gate voltage 
U is used to bias the system such that single-Cooper- 
pair tunneling is resonant across the left J J. Independent 
of the realization of the circuit, the quantized mode in 
a nearby stripline or electromagnetic oscillator has the 
role of the photon cavity and is described as a harmonic 
oscillator {h — 1), -ffphoton = woa^a. 

For the series coupling scheme the interaction Hamil- 
tonian between the atom and the cavity has the form 

i^int = -EjR C0s[2eFi + (t)L+ G{a + a^)] + Fee, (3) 

where Ejn is the Josephson coupling across the right- 
hand JJ and we have used the condition that the phase 
difference across the SSET and the oscillator is fixed by 
the transport voltage V , similar to the approach used in 
Ref. [20I Additionally we have a capacitive coupling term, 
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The important coupling constant is in this case Q = 
{2ec/ELY'\ where ec = e^Ca + Cl + Cr) /2Cl, Cl = 
(Cg + Cl)Ch + C(Cg + Cl + Cr) and El is the mag- 
netic energy of the inductor (the capacitances are defined 
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FIG. 1. (a) The superconducting single-electron transistor 
(SSET) in series with the LC-oscillator. The SSET consists of 
two Josephson junctions (crossed boxes) with a capacitively 
connected gate on the small island in between, (b) Illustration 
of the inductive coupling scheme, where a SQUID forms an 
effective right-hand JJ and is thread by the bias flux (po and 
the oscillator flux iQ{a^ — a), (c) The photon number (n) 
in the oscillator [setup (a)] as a function of the gate charge 
No oc U and the transport voltage V. Two similar resonance 
regions appear corresponding to interchanging the roles of the 
two JJs. A symmetry point of particular interest is circled. At 
this point, photons are created in a cascade of energy decay 
and Cooper-pair tunneling through the right junction (see 
text for more detailed explanation). For our discussion of 
the regime of strongly nonlinear coupling, we restrict ourself 
to this symmetry point. The parameters of the simulation 
are Ejl ^ Ejr ^ 3a;o/10 = 30 ^teV, Cl ^ Cr = 40Cg = 
C/1250 = 0.4 fF (e = 0.08), T = 50 mK, Ru = Rv = 
100 Q., and we have included an additional oscillator Q-factor 
Q = 2 X 10^ 



in Fig. [T]). The charging energy of the SSET island is 
Ec = (C + Cfl,)e^/2C|., the strength of the capacitive 
coupling is characterized by E^c = b'^Cr/C^, and the 
gate charge has the form Na = [V{Cr/{C + Cr)){C + 
Ch{Cg + C)/{Cg + Cl + Cr)) - U{Cg - CICg/{{C + 
Cr){Cg + Cl + C R)))]/2e . Finally, the frequency of the 
oscillator is = \/%tcEL- 

The inductive coupling scheme (Fig. [TJd) results in a 
similar atom-photon coupling. The interaction Hamilto- 
nian is then 



E\„t = EjRsin 
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(5) 
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where we assumed equal Josephson energies of the right- 
hand JJs EjRi = EjR2 = Ejr/2 and an additional 



external flux (po ~ it in the SQUID. Here the coupling 
strength Q depends on the geometry of the system and 
remains small, C/ <C 1, in contrast to the series coupling 
scheme where C/ ~ 1 is possible. For the gate charge we 
obtain Ng = CrV - CgU {Cr = Cri + Cr2) and the 
charging energy becomes Ec — /2{Cl + Cr + Cg)- 

Additionally to the coherent time evolution, incoher- 
ent effects are introduced through voltage fluctuations 
leading io V V + Vf and U U + Uf. It is this 
noise (which is present in any circuit coupled to an elec- 
tromagnetic environment) that creates directionality in 
our lasing cycle, since all other processes are coherent 
and therefore time-reversible. The voltage fluctuations 
are described by the Hamiltonian iJsE and are mediated 
on the SSET by the operators Ny and Nu. The total 
system is then described by a Hamiltonian of the formal 



Ht^ Hq + 2eNvVf + 2eNuUf + Hee ■ 



(6) 



We assume that the voltage fluctuations arc character- 
ized by equilibrium correlations of the form 
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where Ru is a typical line impedance ^ 100 57 with a 
cut off ujc ~ i/RijCu ^ uiQ and that Uf and Vf are 
uncorrelated. For the series LC oscillator one obtains 
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Similar relations apply for transport voltage fluctuations 
Vf and the interaction operator Ny- As one can see for 
series coupling we create additional noise in the oscilla- 
tor as well. We will always include this effect by using 
an oscillator decay rate which will be larger than the de- 
cay induced only by the voltage fluctuations to allow for 
additional internal or external loss. 

For the remainder of the paper, we focus on the se- 
ries coupling scheme (Fig. [T^) as this allows a stronger 
nonlinear coupling between qubit and oscillator, but our 
results apply to both schemes. In Fig. [!}: we plot the 
average photon number as a function of the gate charge 
Ng and transport voltage V for our device, obtained us- 
ing a full density-matrix simulation of the SSET and os- 
cillator. We diagonalize the Hamiltonian ([T]) using the 
Floquet expansion of the eigcnstates. 
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where \n) are the photon number states of the oscillator 
and \N) are the charge states of the island. We then 
treat the voltage fluctuations perturbatively by expand- 
ing the time evolution of the resulting density matrix in 
orders of the coupling to the reservoirs. After tracing out 



the reservoir degrees of freedom we arrive at the Bloch- 
Redfield equation for the reduced density matrix p of the 
quantum system. A detailed discussion of this formahsm 
for the SSET and the methods used to solve the master 
equation can be found in Ref. [2l]. 

Cascade Resonance Several resonances appear at 
different regions in the V — Nq plane shown in Fig. [TJ;. 
These points correspond to a resonance condition be- 
tween levels of the SSET and the oscillator. We focus our 
further discussion on the symmetry point at Nq = 1/2, 
where the charge states |iV = 0) and | = 1) are degener- 
ate. At this symmetry point we are able to qualitatively 
describe features of the system by a simplified Master 
equation. We will now discuss the relevant terms in this 
equation and explain how a lasing cycle can be achieved. 

The eigenstates of the effective atom are given by the 
qubit states 

|t> = ^(|A^-0)-|A^ = l», (10) 

U) = -^(|7V = 0) + |A^ = 1)), 

with the energy splitting /S.E — Ejl. Close to the sym- 
metry point the energy changes with {Nq — 1/2)^ and is 
therefore protected to the first order from low-frequency 
charge noise^^. In the basis of the qubit states we can 
write the total coherent Hamiltonian as 
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+ u!oa^a (11) 
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where we neglected the nonrcsonant capacitivc coupling 
Hcc and at are the Pauli matrices acting on the eigen- 
states (jlOp . Now we perform a transformation U = 
exp [— i2el^a^ai] and, after a rotating-wave approxima- 
tion, this yields 



1 



Hu = -AEaz - ujcsa ' a 
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with an effective oscillator frequency Weg = 2eV — ujo and 
the operators 

s+ =^(n + l|sin[g(at-f a)] \n)\n+\){n\, (13) 



and s+ = sL. At the symmetry point that is circled 
in Fig. [TJ; we have 2eV > loq. This means that in our 
rotating frame each photon that is created effectively de- 
creases the energy of the system. 

The resonance condition is given by loq = 2eV — AE 
(AE = ujoff). In this case there is a coherent transition 
between | 4,)|n) and | ']')\n + 1). To create an overall in- 
crease in the photon number we now need an additional 
incoherent pumping process which is given by the fluctu- 
ations of the electromagnetic environment. In the basis of 
the qubit states (|10p the charge operator is transverse to 
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FIG. 2. The photon distribution p„ as a function of the cou- 
pling strength G, where we keep the linear coupling strength 
EjrQ constant. We see three major types of behavior, each of 
which we show in the inset. A Poisson distribution for small 
coupling (black, J) and for increasing Q a squeezed distribu- 
tion (blue, 77) and multistability (green line, 777). We used 
the parameters: n/y = 0.027 and QEjr/'j = 6.7. 



the eigenstates, N oc <Tx- This means the charge noise al- 
lows a decay from the qubit state | f) to the qubit ground 
state I \.) . This in turn creates a higher population in the 
state I l). Through this mechanism photons are created 
in a cascade of energy decay and absorption of voltage 
quanta. As shown in the schematic in Fig. [TJ;, each time 
a Cooper-pair tunnels through the right junction it cre- 
ates (or absorbs) a photon. The energy to keep the cycle 
going is provided by the voltage source. In the linear 
(small G) limit, this scheme is similar to Rabi-Sideband 
lasing, that is well understood in conventional quantum 
optics^^i^ and has been discussed and demonstrated in 
circuit QEB^. 

At the symmetry point we can write the Master equa- 
tion in a simplified Lindblad form. 



p = -i[Hu,p] + Cchp + CdissP , 
with dissipation in the oscillator 

/^dissP = 2 (2a/0at - [a^a, p]+) , 
and qubit decay 



(14) 
(15) 
(16) 



where [,]-|_ is the anticommutator, and the decay rate is 
given by 

/oo 
dt(t/Kt)C//(0))e'^^* (17) 
-OO 

+ ia|2eA^y|t>P / dtiVfmfme^'^^K 



The overall form of our master equation is now signifi- 
cantly simplified and similar to the master equation of 
a single-atom injection maser. In the rotating-wave ap- 
proximation used to derive the Hamiltonian (|12p we ne- 
glected many matrix elements connecting different pho- 
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ton number states. However, comparison to our full so- 
lution shows that all qualitative features of the system 
are well described by Eq. (fT4|) . 

Anharmonic Coupling The key difference between 
this circuit and existing circuit-QED setups is the anhar- 
monic coupling term, sin [C/(a^' -f a)] , see Eqs. (jlHllSp . 
To understand the effects of this term, we consider three 
different regimes. Expanding the anharmonic term for 
weak coupling Q <^ 1, we obtain the usual constant cou- 
pling term 



sin [g{a^ + a)] = g(at + a) + O ( [g{a + a^)\ 



(18) 



which gives in the lowest order s+ — Qa^ resulting in the 
usual single-qubit lasing results^i^. In this case we find 
the photon number at resonance, Weff = Ai? = Ejl, to 
be2^ 



Wo = 
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for {n)o 3> 1. This expression for the photon number 
in the linear coupling limit will prove to be useful in 
understanding the other regimes that we can reach using 
this device. 

A main advantage of this circuit is that we are not lim- 
ited to the linear coupling regime. As the strength of the 
effective coupling term is roughly sinusoidal with Q and 
photon number, we see a number of new effects. Ironi- 
cally, when operating in this mode, our micromaser be- 
haves in a very similar fashion to a single-atom injection 
maser— ii^ii^, where our oscillatory coupling provides a 
direct analog to the spatially dependent atom-field cou- 
pling in injection masers. In Fig. [5] one can see the prob- 
ability distribution of the photon number states, 



Pn = (n|Tr|/^[p] \n) 



(20) 



in the stationary limit p = 0, as function of the cou- 
pling strength Q. As discussed above, for small coupling 
(marked by /) we simply get the Poisson distribution we 
would expect for a laser—. 

For slightly larger coupling (or photon number) we ob- 
tained a squeezed number distribution within the oscilla- 
tor. This regime (//) is reached when tt/2 < {n)^ < 
TT. If this condition is fulfilled, as t he photon number 
approaches the upper limit (Qy^ (n)o — > tt). the matrix 
element s+ is effectively 'cut off' resulting in the usual 
squeezed state physic a^^'^^ . This can be seen in the asym- 
metric character of the photon number distribution (the 
blue curve) in the inset to Fig. [2] 

Increasing the coupling further {Q\/ {n)o > tt), we 
reach a regime of multistability (HI). In this regime, 
the system can occupy either the original squeezed state, 
restricted by the first zero crossing of sin [t/(a^ -I- a)] , or 
can occupy a new squeezed state associated with the sec- 
ond or subsequent zeros in the sine function. Which of 
these states are occupied depends on which matrix ele- 
ments are small, resulting in 'hot spots' where the system 




FIG. 



3. The amplitude C„ and the phase correlation function 
as a function of time, calculated from a numerical solu- 
tion of the time evolution of the master equation (|14[) . The 
different lines correspond to the different operating regimes 
shown in Fig. [1 (black line I) Q = 0.37, (blue line II) 
Q — 93, (green line ///) Q — 0.968. We used the parameters: 
k/7 = 0.027 and QEjrI^ = 6.7. 



(as a function of ^) suddenly jumps from one squeezed 
state to another, as can be seen in the bright regions 
of Fig. [2l Near the crossover between these states, the 
system displays bistability (and at higher Q multistabil- 
ity) where the oscillator is in a statistical mixture of two 
squeezed states with distinct photon number distribu- 
tions (the green curve in the inset of Fig. [2]). 

Measurement A great advantage of the circuit-QED 
realization of a micromaser is that it is possible to make 
a time resolved measurement of the microwave radiation 
emitted from the LC-oscillator. Measuring the phase and 
the amplitude fluctuations makes it possible to distin- 
guish between the different states of radiation in the cav- 
ity. Amplitude fluctuations are described by the ampli- 
tude correlation function, 



C„(t) = 



(at(t)Q(t)Qt(0)Q(0)) - (flta; 
((ata)2) _ (ata)2 



(21) 



while phase fluctuations are given by the phase correla- 
tion function 



(«Ut)^(0)) 
(ata) 



(22) 



Both correlators have been normalized such that 
Ci(0) = 1 and they approach zero for long times, Ci(t — )■ 
oo) — )• 0. With this normalization, each correlator is fully 
characterized by a single decay rate, Ci(t) oc e"**'*. In 
Fig. [3] we show the amplitude and the phase correlation 
function for the LC oscillator in the various operating 
regimes. One should note that the time- axis is logarith- 
mic since differing regimes produce decay rates that differ 
by several orders of magnitude. 

The two operating regimes // and /// can be clearly 
distinguished in the amplitude correlation function. For a 
squeezed distribution (//), the amplitude is very similar 
at all times and therefore the amplitude decays very fast 
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to its long time limit. If the photon distribution is mul- 
tistabile the situation is very different. In this case the 
field is fluctuating between two favored photon number 
states that are only connected by a small matrix element. 
Therefore it takes a long time until the amplitude corre- 
lator decays to its stationary value. The decay rates of 
the amplitude correlation function are therefore ordered 
such that kIj'^ > > k^^. The decay rate of the coher- 
ent state lies between the two other decay rates as it is 
given approximately by the oscillator decay rate « k. 

To find a clear signature distinguishing the coherent 
state it is necessary to turn to the phase correlation 
function. For the squeezed and the multistable state 
the phase is not well defined and the correlator decays 
quickly, k/J'^^^ ~ k,. In contrast, the coherent state has 
a very long phase correlation time. This is the standard 
result for a laser—, since the phase correlation decay rate 
for a laser is inversely proportional to the average photon 
number, oc l/{n), 

Conclusion We have presented a circuit design for 
creating non-trivial microwave photon distributions us- 



ing a SSET strongly coupled to a strip-line or LC res- 
onator. Using a combination of Cooper-pair tunneling 
and coupling to voltage fluctuation noise, a lasing cycle 
can be established. As this can be achieved while op- 
erating at a charge degeneracy point, the system is also 
protected against low-frequency charge noise. 

Realizing the strongly nonlinear coupling regime re- 
quires an oscillator with a large ratio of charging energy 
to inductive energy, Q ~ {2ec/ELy^^. Such an oscil- 
lator has been demonstrated in Ref. [2^. A supercon- 
ducting quantum interference device array was used as 
an effective resonator, which has the additional advan- 
tage of allowing for a tunable Q, such that all regimes 
discussed in this paper can be accessed in a single ex- 
periment. Since the overall photon number can remain 
small the anharmonicity of the SQUID array should not 
play significant role. Another possible realization is to 
use a tunable right-hand (or left-hand) Josephson junc- 
tion. This would allow to change the photon number of 
the linear limit (n)o. An increase of the photon number 
would allow the measurement of all three regimes as well. 
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